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Abstract 

Noise and decoherence are major obstacles to the implementation of Josephson junction qubits in 
quantum computing. Recent experiments suggest that two level systems (TLS) in the oxide tunnel 
barrier are a source of decoherence. We explore two decoherence mechanisms in which these two 
level systems lead to the decay of Rabi oscillations that result when Josephson junction qubits 
are subjected to strong microwave driving. (A) We consider a Josephson qubit coupled resonantly 
to a two level system, i.e., the qubit and TLS have equal energy splittings. As a result of this 
resonant interaction, the occupation probability of the excited state of the qubit exhibits beating. 
Decoherence of the qubit results when the two level system decays from its excited state by emitting 
a phonon. (B) Fluctuations of the two level systems in the oxide barrier produce fluctuations and 
1/f noise in the Josephson junction critical current I Q . This in turn leads to fluctuations in the 
qubit energy splitting that degrades the qubit coherence. We compare our results with experiments 
on Josephson junction phase qubits. 

PACS numbers: 03.65.Yz, 03.67.Lx, 85.25.Cp 
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I. INTRODUCTION 



The Josephson junction qubit is a leading candidate as a basic component of a quantum 
computer. A significant advantage of this approach is scalability, as these qubits may be 
readily fabricated in large numbers using integrated-circuit technology. Recent experiments 
on Josephson qubits have successfully shown that they possess quantum coherent properties 



QflflflflflflQ. 



However, a major obstacle to the realization of quantum computers 
with Josephson junction qubits is decoherence and the dominant noise source has not yet 
been identified. Because the measured decoherence times are substantially shorter than 
what is needed for a quantum computer, there has been ongoing research to understand 
decoherence mechanisms in Josephson qubits. Let us briefly review the theoretical work 
that has been done. Martinis et al. jj| investigated decoherence in Josephson phase qubits 
due to current noise primarily from external sources. Paladino et al. |t| analyzed decoherence 



in charge qubits due to background charge fluctuations. Van Harlingen et al. [10|, studied 
how low frequency 1/ f critical current fluctuations lead to decoherence in various types of 

n 

Josephson junctions. (/ denotes frequency.) Smirnov [12| used perturbation theory to study 
the decay of Rabi oscillations due to thermal fluctuations of a heat bath weakly coupled to 
the qubit. Even though fluctuating two level systems in oxide tunnel barriers have long 
been known to be a major intrinsic noise source in Josephson junctions jlUl. 11^. their role 
in qubit decoherence has not been investigated theoretically. In this paper we will focus on 
the intrinsic microscopic mechanisms whereby two level systems produce decoherence in a 
Josephson qubit. 

A recent experiment indicates that two level systems can couple to a Josephson phase 
qubit In these experiments the lowest excitation frequency of a qubit is measured as 
a function of the bias current that determines the depth of the potential energy double 
well of the qubit. For most values of the current bias, a single excitation frequency Uiq is 
observed and it decreases with increasing bias current. Occasionally the experiments find 
spurious resonances characterized by two closely spaced excitation frequencies at a given 
bias current. The size of the gap between the two excitation frequencies is on the order of 



25 MHz. Simmonds et al. 



15| have argued that this splitting is evidence that the qubit 



is coupled to a two level system (TLS) with an energy splitting very close or equal to that 
of the qubit. We will refer to this as a resonant interaction. The microscopic nature of 
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these two level systems is unclear. It may be due to the motion of oxygen atoms in the 



oxide tunnel barrier of the Josephson junction [16. 1171 ] . An oxygen atom (or vacancy) could 
sit in a double well potential and tunnel between two positions. Or it may be due to a 
quasiparticle hopping between two positions in the oxide barrier. Or an electron trap in 
the oxide barrier could fluctuate between being occupied and empty 3, 14, [lfl 18]. Or a 
trapped flux quantum could be tunneling back and forth between two positions in the oxide 
barrier. Regardless of the microscopic nature of the TLS, we can use the fact that it is a 
two level system to understand how such a defect can couple to the qubit. Note that the 
qubit energy splitting is a function of the critical current I Q . Fluctuations of the TLS lead 
to fluctuations of the tunneling matrix element T through the oxide barrier. This in turn 
leads to fluctuations in the critical current I Q since I D ~ |T| 2 . If the two states of the TLS 
correspond to two different values, J 0j i and J 0j 2, of the critical current, the coupling between 
the qubit and the TLS is proportional to the difference |J i — J 0>2 |. 

Experiments often probe qubits using Rabi oscillations. Let us take a moment to review 
Rabi oscillations Q]. If the qubit is initially in its ground state, resonant microwaves with a 
frequency that matches the qubit energy splitting (c^io) will initially increase the probability 
amplitude of finding the qubit in its excited state (|1)). However, as time goes on, at some 
point the qubit is completely in its excited state, and the electromagnetic wave goes on to de- 
excite the qubit through stimulated emission. Thus the system will be coherently oscillating 
between the two energy eigenstates with a Rabi frequency /r. The frequency fn of the Rabi 
oscillations increases linearly with the amplitude of the driving electric field. Rabi oscillations 
lave been seen in the occupation probability Pi of the excited state of a Josephson qubit 
r|. This demonstration of quantum coherence is a preliminary requirement for quantum 
computing but most of the reported Rabi oscillations in Josephson qubits have a rather 
small amplitude (less than 50%) and a short coherence time (less than one microsecond). 
Experiments have found that the presence of a resonant interaction between the qubit and 
a two level system substantially reduces or even eliminates Rabi oscillations jisl ]. 

In this paper we theoretically model a qubit coupled to a two level system and study 
the effect of this coupling on Rabi oscillations. We study the quantum dynamics of a 
Josephson qubit by numerically integrating the time-dependent Schrodinger equation. Our 
method is not limited to weak coupling to the noise sources. While our analysis applies 
to any Josephson qubit, for illustration we consider the Josephson phase qubit that was 
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studied in recent experiments |15l |. We explore two decoherence mechanisms where two level 
systems lead to the decay of Rabi oscillations that result when Josephson junction qubits 
are subjected to strong microwave driving. We first consider the resonant case. In the 
resonant regime, the energy splitting of the two level system and the qubit are matched. 
As a result the occupation probability of the qubit 's excited state exhibits beating. This 
has been termed a qubit duet [?], 3|. Our calculations show that decoherence of the qubit 
results when the two level system decays from its excited state by emitting a phonon. In 
section II, we numerically calculate the qubit occupation probability as a function of time 
and compare it with experiment. 

The other case involves low frequency fluctuations of the qubit energy splitting. Even 
though the qubit energy splitting u w / '27i is on the order 10 GHz, low frequency fluctuations 
of the critical current lead to low frequency fluctuations of the qubit energy splitting. We 
hypothesize that this noise comes from slow fluctuations of two level systems in the oxide 
tunnel barrier. If we assume that there are a number of two level systems in the barrier, and 
if these two level systems have a broad distribution of decay rates, then they will produce 
1/ / noise j^fll. I21I] that leads to decoherence. In section III we show that these fluctuations 
in the qubit energy splitting lead to decay of the Rabi oscillations. We consider three cases. 
In the first case the qubit is coupled to a single slowly fluctuating two level system. In the 
second case the qubit is coupled resonantly to a two level system and, at the same time, is 
subjected to slow fluctuations in the critical current of a two level system. In the third case 
we consider a qubit with energy splitting fluctuations that have a 1/ f noise spectrum. 



II. QUBIT-TLS RESONANCE 



In our model we assume that there are two level systems in the oxide barrier of the 



22 



was 



n 



Josephson junction. The standard model of noninteracting two level systems 
introduced by Anderson, Halperin, and Varma [^J, and independently by W. A. Phillips |2ifl 
in 1972. The standard Hamiltonian for a two level system is 

^ \ V.I 

(1) 
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Here we are using the left well - right well basis where \L) (\R)) is the left (right) well state. 
A is the asymmetry energy, i.e., A is the energy difference between the right well and the 



left well. We can diagonalize the Hamiltonian to get the energy eigenvalues that are given 
by ±etls/2 where 



e TLS = VA 2 + A 2 . (2) 

Simmonds et al. have found experimental evidence for resonant interactions between a 
phase qubit and a two level system Q|. This resonant interaction occurs when the energy 
splitting eTLS of the TLS matches that of the qubit, i.e., €tls — ^io- Simmonds et al. 
constructed a phenomenological model to account for their experimental findings |15| . In 
this section we investigate Rabi oscillations in the presence of a qubit-TLS resonance. We 
incorporate decoherence into the model of Simmonds et al. to show that Rabi oscillations 
can decohere due to decay of the excited state of the two level system via phonon emission. 
We briefly describe the model of Simmonds et al. in the following. 

The Hamiltonian of a Josephson phase qubit (which is essentially a current-biased Joseph- 



son junction) is 2f 



27, 



where $ = h/2e is a superconducting flux quantum. C is the capacitance of the junction, 
and hias is the bias current. The operators Q and 5 correspond to the charge and phase 
difference across the Josephson junction respectively. 

Next we assume that there is a two level system in the barrier of the Josephson junction. 
The two states of the TLS correspond to two different values of the Josephson junction 
critical current I a which is proportional to the square of the tunneling matrix element. 
When the TLS is in state \R) (state |L)), the junction critical current is I r (Iql)- The 
qubit couples to the TLS because the qubit 's energy splitting huio is a function of I Q . The 
expression for uiq can be derived from the Hamiltonian in Eq. (J3J) using the expression for 
the resonant frequency of an LC circuit ujio « 1/y/LjC where Lj = Q /2irI cos5 is the 
Josephson inductance 2^ |: 

1/4 
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(4) 



Here we used the fact that / = I a smS implies that cos<5 = J 1 — (I/I ) 2 . Typically hias is 
slightly less than I a . The qubit couples to the TLS because the Josephson junction critical 
current I Q is modified by the TLS. Therefore, the interaction Hamiltonian is 

H qb -TLs = -^^cos5®\R) (R\ - %^cos(5® \L) (L\ . (5) 



We can transform to the eigenbasis of the TLS. Provided that the TLS is symmetric, its 
ground state is \g) = (\R) + \L))/y/2 and its excited state is |e) = (\R) — \L))/y/2. The 
ground state of the qubit is |0) and the excited state of the qubit is |1). We can rewrite the 
operator cos 5 in terms of its matrix elements in the qubit basis by noting that the phase 
qubit is typically biased close to 5 = tt/2 where the Josephson current is maximized. So 
cos(7r/2 — 5') = sin(<5') ~ 5' where i5'C 1. 5' can be represented as a sum of a creation and 
an annihilation operator in much the same way as the position coordinate of a harmonic 
oscillator. Using these facts, we find that Eq. (J3J) becomes |3| 
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Hqb-TLS = -yy 2^ GM(M + |l,e)<0,<?| + |l,</)<0,e| + |0,e)(l,^|) 



(6) 



where uj\q is the energy difference in the qubit levels, and 5I = Iqr — Iql is the fluctuation 



can be 



amplitude in I Q produced by the TLS. Since the values of (H qb _TLs), C and ui 
determined experimentally, one can estimate 5I /I to be approximately 6 x 10~ 5 15| when 
the Josephson junction is in the zero- voltage state. However this value of 5I /I differs from 
previous measurements made on Josephson junctions in the finite voltage state that found 
low frequency (1 kHz or less) fluctuations with 8I /I ~ 2 x 10~ 6 HQ. 

Next we study the quantum dynamics of the coupled qubit-TLS system when it is sub- 
jected to microwave driving at a frequency ujiq equal to the qubit energy splitting. We 
will assume that both the qubit and the TLS can couple to microwaves. The resulting 
Hamiltonian matrix of the qubit-TLS model is, 
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where the basis states are |0, g), |0,e), \l,g), and |l,e) respectively, r] = (5I o /2)yh/(2ui C) 
is the coupling between the qubit and the TLS; g qb is the coupling between the qubit and 
the microwaves; and gris is the coupling between the TLS and the microwaves. In the 
strong driving regime, the qubit-microwave coupling is larger than the qubit-TLS coupling. 
We consider the case of strong driving throughout the paper because it is the experimental 
condition in Ref . . Without microwave driving, the Hamiltonian matrix can be decoupled 
into two 2x2 matrices. We calculate the time evolution of the qubit-TLS wave function by 



integrating the time-dependent Schrodinger equation using the Runga-Kutta method j^. 
The initial wave function \ip) is the ground state of Eq. ((7|). When huio and Etls are much 
greater than 77, the ground state wave function is approximately given by |0, g). 

We first consider the case of strong driving with gxis — and with the TLS in resonance 
with the qubit, i.e. Etls — ^10 ■ If there is no coupling between the qubit and the TLS, 
then the four states of the system are the ground state |0, g), the highest energy state |1, e), 
and two degenerate states in the middle \l,g) and |0,e). If the qubit and the TLS are 
coupled with coupling strength 77, the degeneracy is split by an energy 2r). Figure [T] shows 
the coherent oscillations of the resonant qubit-TLS system. We define a projection operator 
Pi = \l,g)(l,g\ + 1 1, e) { 1, e | so that (Pi) corresponds to the occupation probability of the 
qubit to be in state |1) as in the phase-qubit experiment. Instead of being sinusoidal like 
typical Rabi oscillations (the dotted curve), the occupation probability Pi exhibits beating 
(Fig. la) because the two entangled states that are linear combinations of \l,g) and |0,e) 
have a small energy splitting 2i], and this small splitting is the beat frequency. Without any 
source of decoherence, the resonant beating will not decay. Thus far the beating phenomenon 
has not yet been experimentally verified. The lack of experimental observation of beating 
implies that the TLS or qubit decoheres in less time than the period ~ I/77. Note in Fig. ^1 
that the second beat is out of phase when compared with the usual Rabi oscillations. The 
occupation probabilities in the individual states are plotted in Fig. lb-le. (To the best of our 
knowledge, these quantities are not measurable.) We find that there is a very low occupation 
probability in the states |0,e) and |l,e) during the first Rabi cycle because the TLS is not 
directly coupled to the microwaves, and thus the TLS tends to be in its ground state \g). 
In the limit 77 — > 0, the system oscillates coherently between \0,g) and \l,g). Occupying 
the two states |0,e) and |l,e) occurs only via the qubit-TLS resonance coupling 77. From 
Fermi's golden rule, the average transition rate to |0,e) from \l,g) is 2irr] 2 /h which is much 
slower than the initial transition rate from \0,g) to \l,g). 

Now we consider the decay of excited TLS via phonons as a source of decoherence for the 
qubit-TLS system. Two level systems couple to the strain field and are able to decay from 
their excited state by emitting a phonon with an energy equal to the TLS energy splitting 



The rate for an excited TLS to emit a phonon and return to its ground state is given by 
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where 7 is the deformation potential, p is the mass density, c\ (c t ) is the longitudinal (trans- 
verse) speed of sound, and (3 is the inverse temperature. From Eq. (jHJ), the relaxation time 
T p h is estimated to be in the range from 10 to 100 ns. (We find r p h ~ 80 ns when using the 
following values that are appropriate for a symmetric TLS in SiC>2 with a tunnel splitting 
that matches the energy splitting of a Josephson junction qubit: 7 = 1.0 eV, p = 2.2 g/cm 3 , 
q = 5.8 x 10 5 cm/s, c t = 3.8 x 10 5 cm/s, A = 0, A G = 0.5K, and T = 25 mK.) 

The decay of the excited two level system leads to decoherence of the Rabi oscillations. 
We can incorporate this relaxation rate of the excited TLS into our calculations of the Rabi 
oscillations of the qubit-TLS system by using the Monte Carlo wave function method j^j]. 
In the original application |3j| a two-level atom driven by a laser field can decay by emitting 
a photon. We can easily generalize the algorithm for our somewhat more complicated case 
in which the qubit-TLS system decays either from |0, e) to |0, g) or from |1, e) to |1, g). The 
algorithm goes as follows, (a) Numerically propagate the wave function from time tj to 
ti + At as if there is no energy decay in the TLS. Note that the probability for decaying 
from |0,e) to |0,g) during the time interval At is P(o, e ) x 1 — ex P (777 ) ' wnere ^(o.e) is 
the probability that state |0,e) is occupied, (b) Generate a uniformly distributed random 
number G [0, 1]. If r.; < P(o,e) x 1 — exp (7—^) , then the system decays from |0,e) to 



T ph ) _ 

|0,g), and we represent this by resetting the wave function ip(U + At) to |0, g). Then we 
repeat steps (a) and (b). Similarly, if 1 — < P(i, e ) x 1 — exp (77^) , then the system 
decays from |l,e) to \l,g), and we represent this by resetting the wave function ip{ti + At) 
to |1,(7). Then we repeat steps (a)-(b). Notice that 1 — exp (— — ) is very small since the 



Tph /. 

time step At is small. So there is no chance that both decays could happen in the same 
time step. If neither of the above criteria are satisfied, then the qubit does not decay. Then 
we repeat steps (a)-(b) and keep propagating the wave function until the desired finishing 
time. 

We have used our algorithm to study the effect of TLS energy decay via phonon emission 
on Rabi oscillations. The results are shown in Fig. |21 The dotted lines show the beating that 
occurs when the qubit and TLS are in resonance with no TLS decay. The solid lines shows 
the rapid damping of the Rabi oscillations and the dephasing that occurs when the TLS can 
decay via phonon emission. In Fig. |2k there is no direct coupling between the microwaves 
and the TLS, whereas in Fig. EJd, the microwaves are directly coupled both to the qubit 
and to the TLS. We see that in the latter case beating is damped out more quickly. The 
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Rabi decay time trcm can be denned as the time for the envelope of the Rabi oscillations to 
decay by 1/e of their original amplitude. Fig. El shows that trum is longer than r p h simply 
because the excited state of the TLS is not always occupied and available for decay. Fig. ^ 
shows that states |0, e) and |1, e) are essentially unoccupied at early times and that the TLS 
energy decay can take place only when the excited state of the TLS is sufficiently populated. 
Therefore in Fig. the solid and dotted lines are in phase for the first few Rabi cycles. 

Figure Eh> shows the effect of coupling between the microwaves and the TLS. It has 
been shown experimentally that a TLS will couple to microwaves if it has an intrinsic or 
induced dipole moment 0,0]. We can estimate the value of the matrix element Qtls to be 
approximately pE where p is the magnitude of the dipole moment and E is the magnitude 
of the electric field produced by the incident microwaves. We can estimate p by assuming 
that the TLS is a charged particle (with charge \q\ = e) hopping between two sites separated 
by a couple of angstroms. We can estimate E ~ V/d where d ~ 20 A is the thickness of 
the oxide barrier and V is the voltage produced across the junction by microwaves. We can 
calculate V using the Josephson equations with the current being the microwave current I^ w 
across the junction. Using the experimental values from Ref. jl5l | . we can estimate I^ w ~ 1 
pA. Therefore, we have Qtls ~ 0.017iu;io, which is comparable to g q t,. Comparing the two 
panels, we find that the TLS-microwave coupling greatly enhances the resonant decoherence 
mechanism because the transition rates from the ground state to the states |0,e) and |1, e) 
increase as gris increases. We conclude that the energy decay of the TLS mainly causes 
the Rabi oscillations to decay, and adding the TLS-microwave coupling further degrades 
the qubit coherence. However, the phase qubit experiment indicates that the resonant 
interaction affects both the Rabi amplitude and decay time 3^. The large amplitude of the 
Rabi oscillations that we see at short times is not seen experimentally. Experimentally, a 
qubit in resonance with a TLS has very small amplitude Rabi oscillations at all times. This 
implies that our calculations do not include all the sources of decoherence responsible for 
the experimental observations, such as resonance with multiple fluctuators and interactions 
between the fluctuators. 

Next we propose an experiment to test if the TLS in the Josephson qubit couples to 
microwaves or not. Start in the ground state, then send in a resonant microwave 7r-pulse 
and stop pumping. (A 7r-pulse lasts for half of a Rabi cycle.) When Qtls ^ 9gb, the 
qubit-TLS system mainly occupies the state |1, g) right after the 7r pulse. State \l,g) is a 
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superposition of two eigenstates = (|0, e) + \l, g))/\/2 and \ij;' 2 ) = (|0, e) — 11, g))/\/2) of 
the qubit-TLS system. Thus the wave function oscillates between the states 1 0, e) and |1, g) 
after the -zr-pulse. Provided that the relaxation time of the TLS is longer than 1/rj, coherent 
oscillations in the occupation probability Pi can be observed with a period of 2rj |35[. If 
there is no energy decay of the excited TLS, there is no mechanism for decoherence and the 
oscillation amplitude is one. When Qtls is comparable to g g b, both states \l,g) and |l,e) 
are partially occupied right after the 7r-pulse. State |l,e) is essentially an stationary state. 
Thus partially occupying |1, e) reduces the Rabi oscillation amplitude of Pi and merely gives 
a constant contribution to Pj_. Figure El shows the dynamics of the qubit-TLS system during 
and after a microwave 7r-pulse. As we expect, Fig. Ek shows that the oscillation amplitude 
decreases as gxLS increases. In Fig. Et the excited TLS has no means of decay (r ph = oo), 
so the coherent oscillations do not decay. In Fig. Eb and Fig. Et, we set r p h = 40 ns. As 
a result, the oscillatory P 1 is attenuated. We include the effect of the energy decay of the 
TLS by using the Monte Carlo wave function method described previously. These model 
calculations suggest that one can estimate gris by measuring the oscillation amplitude. An 
alternative way to measure gris is to compare the Rabi oscillations after 7r and 2>n pulses. 
If grLS = 0, right after a ir or 3ii pulse, the system primarily occupies the state |1, g) with 
P(i, g ) ~ 1- Therefore the coherent oscillations after the ir and 3n pulses ought to be similar. 
However, if gris is nonzero, the longer pulse will pump more weight into state |0,e). So 
the amplitude of the oscillations of Pi right after a 7r-pulse will be greater than immediately 
after a 37r-pulse. With a 37r-pulse, the oscillations in Pi are more damped, as shown in 
Fig. Efc. 

III. THE LOW-FREQUENCY TWO-LEVEL FLUCTUATOR 

As described in Sec. I, a two level system trapped inside a Josephson junction barrier 
can produce noise in the critical current I Q by varying the height of the tunneling potential 
barrier. As a result, the qubit energy levels fluctuate and this leads to decoherence. We 
now study the decoherence produced by a low-frequency TLS. As we mentioned earlier, the 
dependence of u 1Q on I a is given by Eq. HI We note that u w is modulated as I a varies. As the 
TLS modulates I Q , the phase of the qubit is randomized and coherent temporal oscillations 
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are destroyed. Using the fact that hias is slightly smaller than I Q in Eq. (JJJ), one finds 

8uj w 8I a 

(u W ) 4 ((J ) - I b ias) 



(9) 



In Ref. 15], it was found that 81 /(I ) ~ 6 x 10 5 . In addition, the phase qubit is typically 



operated at a bias current such that ((J ) — hias)/ (h) — 0.0025. Substituting these numbers 
into Eq. (JHJ), one can estimate that the amplitude of the level fluctuations 8uiq/{uj\q) is 
approximately 0.006. 

The Rabi oscillations are calculated by integrating Schrodinger's equation in the pres- 
ence of noise from a single low- frequency fluctuator, i.e., u>i t TLS < 1. txLS is the char- 
acteristic time of the random switching of the TLS. We simulate the noise using the 
Monte Carlo method. At each time ti, a random number £ [0, 1] is generated. If 
r\ < ((time step) /txLs) < 1, then the two level system switches wells. We can imagine 
that these are thermally activated transitions between the two wells of a symmetric double 
well potential. We assume here that the dwell times £tls are the same in the two wells. 
Each time the TLS switches wells, the critical current, and hence U\q switch between two 
values. The Rabi oscillations of the qubit are calculated using the Hamiltonian 

/ 
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g qb sin(uj w t) 
g qb sm{u 1Q t) h[u w + 8u w (t)} 



(10) 



where the qubit energy levels (|0) and |1)) are the basis states and the noise is produced by 
a single TLS. Our calculations are oriented to the experimental conditions and the results 
are shown in Fig. 0] In Fig. EJi-c the characteristic fluctuation rate t^ l LS = 0.6 GHz. Panel 
shows that the qubit essentially stays coherent when the level fluctuations are small 
(8(Jio/(Jiq — 0.001). Panel H^i shows that when the level fluctuations increase to 0.006, the 
Rabi oscillations decay within 100 ns. The Rabi relaxation time also depends on the Rabi 
frequency as panel 0£ shows. The faster the Rabi oscillations, the longer they last. This is 
because the low-frequency noise is essentially constant over several rapid Rabi oscillations 



12j. Alternatively, one can explain it by the noise power spectrum Si(f). Since the noise 
from a single TLS is a random process characterize d by a single characteristic time scale 



^tls-i it has a Lorentzian power spectrum 
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We do not expect Rabi oscillations to be sensitive to noise at frequencies much greater than 
the frequency of the Rabi oscillations because the higher the frequency /, the smaller the 
noise power and because the Rabi oscillations will tend to average over the noise. Rabi 
dynamics are sensitive to the noise at frequencies comparable to the Rabi frequency. In 
addition, the characteristic fluctuation rate plays an important role in the rate of relaxation 
of the Rabi oscillations. It has been shown that t^is can be thermally activated j3] for 
TLS in a metal-insulator-metal tunnel junction. If the thermally activated behavior applies 
here, the decoherence time r Rabi should decrease as temperature increases. In Fig. 0Jl, the 
characteristic fluctuation rate has been lowered to 0.06 GHz (which is much lower than 
uio/2-k m 10 GHz). The noise still causes qubit decoherence but affects the qubit less than 
in Fig. Fig. E] shows that the noise primarily affects the Rabi amplitude rather than the 
phase. 

Experimentally, the two TLS decoherence mechanisms (resonant interaction and low- 
frequency level fluctuations) can both be active at the same time. We have calculated the 
Rabi oscillations in the presence of both of these decoherence sources by using the qubit-TLS 
Hamiltonian in eq. (jZJ) with a fluctuating uio{t) that is generated in the same way and with 
the same amplitude as in Figure^. We show the result in Fig. El By comparing Fig. El with 
Fig. 2b, we note that adding level fluctuations reduces the Rabi amplitude and renormalizes 
the Rabi frequency. The result in Fig. El is closer to what is seen experimentally . 

In a charge qubit, the qubit energy splitting is M /2e when biased at the degeneracy 
point [2] . So for a charge qubit critical current fluctuations are related to fluctuations in the 
qubit energy splitting by 

5u w = 5Ip 

(«10> (la) ' 1 ] 

Comparing this to the analogous expression for a phase qubit (Eq. (jHJ)), it is obvious that for 
the same amount of noise in the critical current 5I /(I ), the level fluctuations 5uj w /u)io of 
a charge qubit seem to be less than those of a phase qubit. This is because the phase qubit 
must be biased with a current slightly less than (I Q ), which magnifies the critical current 
noise typically by a factor of 100. However, the area of the Josephson junctions in charge 
qubits is much smaller than in phase qubits. As a result , (I n ) is smaller in the charge qubit. 
(In Ref. |2j, the charge qubit has I a ~ 30 nA. In Ref. the phase qubit has I a ~ 11.6 
yiiA.) Thus for a fixed value of SI , the charge qubit has a much larger relative fluctuation 
SI /(I ). Therefore we believe that level fluctuations remain a concern for charge qubits. 
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1/f noise in the critical current of Josephson junctions has been observed experimentally 
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44j. The microscopic source of this noise is unknown, though it has been 
suggested that the tunneling of atoms or ions is involved [3]. When an ensemble of two 
level systems with a distribution of relaxation times is considered, we can obtain the resulting 
noise power spectrum by averaging the Lorentzian power spectra of individual TLS over the 



distribution of relaxation times |20|. For example we can replace tris i n the Lorentzian in 
Eq. (JTT|) with the TLS relaxation rate ttls given by Eq. (jHJ), and average over the TLS 



parameters A and A 



2l|. The result is a 1/f noise spectrum in the critical current, and 



hence in the qubit level fluctuations 8w\q. So we have considered the effect of 1/f noise in a 
Josephson junction on qubit decoherence. We do not average over Lorentzians to obtain our 
noise spectrum. Rather we numerically generate a time series 5u\q with a 1/f noise power 
spectrum. We use this 8u\q in Eq. (fTUj) and calculate the Rabi oscillations with the same 
numerical approach as in Fig. HJ The results are shown in Fig. |B] When the noise power 
is low, the Rabi amplitude is reduced but it remains in phase with the unperturbed Rabi 
oscillations (dotted line). However, as Fig. shows, the coherent nature will eventually be 
destroyed as the noise power increases. The decaying Rabi oscillations clearly demonstrate 
that 1/f noise is able to adversely affect the coherence of the qubit. 



IV. CONCLUSION 



We have shown that the coupling between a two-level fluctuator and a Josephson qubit 
has a large effect on the decay and decoherence of the Rabi oscillations. We have studied 
the quantum dynamics of a Josephson qubit subjected to microwave driving by numerically 
integrating the time-dependent Schrodinger equation. We focussed on two decoherence 
mechanisms of a Josephson qubit caused by coupling to two level systems. (A) In the 
resonant regime we considered a qubit in resonance with a high-frequency TLS, i.e., a TLS 
with a large tunneling splitting. Without any source of energy decay, the qubit occupation 
probability Pi of the excited state exhibits beating rather than being sinusoidal as in the 
usual Rabi oscillations. Including the energy decay of the excited TLS via phonon emission 
results in the decoherence of the coupled qubit-TLS system. This decoherence mechanism 
primarily causes a characteristic relaxation time of the Rabi oscillations. The Rabi dynamics 
at short times is not affected. Furthermore, coupling between the TLS and the microwave 
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driving is found to further degrade qubit coherence. (B) The other regime involved qubit 
level fluctuations. Low-frequency TLS are treated as noise sources because they randomly 
modulate the junction critical current I Q . Fluctuations in I a modulate the qubit energy 
splitting, thus randomizing the phase of the qubit, leading to decoherence. Based on noise 
measurements in prior experiments, our model calculations suggest that noise from a single 
TLS can cause Rabi oscillations to decay within 100 ns. When the qubit is coupled to 
a single slow fluctuator, we have shown that the Rabi decay time depends on the noise 
amplitude 5cu w , the characteristic fluctuation rate, and the Rabi frequency. When the qubit 
level fluctuations have a 1/ f noise spectrum, the Rabi oscillation degradation increases with 
increasing noise power. 
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time (ns) 

FIG. 1: Rabi oscillations of a resonantly coupled qubit-TLS system with etls = fofio- There is 
no mechanism for energy decay. Occupation probabilities of various states are plotted as functions 
of time, (a) P\ is the occupation probability in the qubit state |1); (b) -P(o,g) is the occupation 
probability in the state |0,g); (c) -P(o, e ) is the occupation probability of the state |0, e) ; (d) P(\ g ) is 
the occupation probability of the state \l,g); and (e) P(i e ) is the occupation probability of the state 
1 1 , e) . Notice the beating with frequency 2r). Throughout the paper, loiq/2tt = 10 GHz. Parameters 
are chosen mainly according to the experiment in Ref. Q]: r]/fruio = 0.0005, g q b/tiwio = 0.01, 
and gTLS = 0. The dotted line in panel (a) shows the usual Rabi oscillations without resonant 
interaction, i.e. rj = 0. 



18 




time (ns) 

FIG. 2: Relaxation of Rabi oscillations due to energy decay of the TLS via phonon emission. Solid 
lines include the effect of the energy decay of the TLS with r p h = 10 ns and the dotted lines have 
T p h = oo. These results are averaged over 100 runs in this figure and in all of the following figures, 
(a) No coupling between the TLS and the microwaves, (b) Direct coupling between the TLS and 
the microwaves with qtlf/^w = 0.008. The rest of the parameters are the same as in Fig.l. We 
note that the coupling between the TLS and microwaves degrades the qubit coherence. 
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time (ns) 

FIG. 3: The qubit-TLS system starts in its ground state at t = 0. A microwave tt or 3ir pulse 
(from to 5 ns) puts the qubit-TLS system in the qubit excited state |1) that is a superposition 
of the two entangled states \ip[) = (|0, e) + \l,g))/y/2 and \ip' 2 ) = (|0,e) - \l,g))/y/2. After the 
microwaves are turned off, the occupation probability starts oscillating coherently. Values of gTLS 
indicated in the figure are normalized by huio. The rest of the parameters are the same as in Fig. 
1. (a) No energy decay of the excited TLS, i.e., r p h = oo. Coherent oscillations with various values 
°f 9tls- (b) Oscillations following a -zr-pulse with r p h = 40 ns and various values of gTLS- ( c ) 
Oscillations following a 7r-pulse and a 37r-pulse with t p h = 40 ns and gTLS = 0.004. The dip in the 
dot-dash line is one and a half Rabi cycles. 
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FIG. 4: Solid lines show the Rabi oscillation decay due to qubit level fluctuations caused by a 
single fluctuating two level system trapped inside the insulating tunnel barrier. The TLS produces 
random telegraph noise in I Q that modulates the qubit energy level splitting (a) The level 

fluctuation 5uj\q/{ujiq) = 0.001. The characteristic fluctuation rate t^,\ s = 0.6 GHz. The Rabi 
frequency fn = 0.1 GHz. The dotted lines show the usual Rabi oscillations without any noise 
source, (b) <5a;io/(wio) = 0.006, t^ l LS = 0.6 GHz, and fn = 0.1 GHz. The dotted lines show the 
usual Rabi oscillations without any noise source, (c) 8uiq/{uj\q) = 0.006, t^\ s = 0.6 GHz, and 
f R = 0.5 GHz. (d) <Vl /(wio) = 0.006, t^ LS = 0.06 GHz, and f R = 0.5 GHz. Note that the scales 
of the horizontal axes in (a)-(c) are the same. They are different from that in (d). 
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FIG. 5: Solid line represents Rabi oscillations in the presence of both TLS decoherence mechanisms: 
resonant interaction between the TLS and the qubit, and low frequency qubit energy level fluctu- 
ations caused by a single fluctuating TLS. The TLS couples to microwaves (gTLS /(^io) = 0.008) 
and the energy decay time for the TLS is r p h = 10 ns, the same as in Fig. 2b. The size of the 
qubit level fluctuations is the same as in Fig. \$jp. The dotted line shows the unperturbed Rabi 
oscillations. 
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FIG. 6: (a)-(b): Rabi oscillations in the presence of 1/f noise in the qubit energy level splitting. 
Dotted curves show the Rabi oscillations without the influence of noise. Panel (c) shows the two 
noise power spectra S(f) = |&t>io(/)/(wio)| 2 of the fluctuations in u>iq that were used to produce 
the solid curves in panels (a) and (b). Rabi frequency /r = 0.1 GHz. 
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